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STABILITY UNDER GALERKIN TRUNCATION OF A-STABLE 
RUNGE-KUTTA DISCRETIZATIONS IN TIME 

MARCEL OLIVER AND CLAUDIA WULFF 



Abstract. We consider semilinear evolution equations for which the Hnear 
part is normal and generates a strongly continuous semigroup and the nonlin- 
ear part is sufficiently smooth on a scale of Hilbert spaces. We approximate 
their semiflow by an implicit, A-stable Runge-Kutta discretization in time and 
l/^ ■ a spectral Galerkin truncation in space. We show regularity of the Galerkin- 

truncated semiflow and its time-discretization on open sets of initial values 

with bounds that are uniform in the spatial resolution and the initial value. 

We also prove convergence of the space-time discretization without any con- 

^r ' dition that couples the time step to the spatial resolution. Then we estimate 

►^^ ' the Galerkin truncation error for the semiflow of the evolution equation, its 

^^ , Runge-Kutta discretization, and their respective derivatives, showing how the 

r^ ■ order of the Galerkin truncation error depends on the smoothness of the initial 

data. Our results apply, in particular, to the semilinear wave equation and to 
the nonlinear Schrodinger equation. 



(N. 

K* " Contents 

<N 

1. Introduction 1 

C^^ ■ 2. Semilinear evolution equations under Galerkin truncation 4 

~~J, . 2.1. General setting 4 

C^ ■ 2.2. Spectral Galerkin truncation and convergence 6 

^3) , 2.3. Regularity of Galerkin truncated semiflow 8 

2.4. Accuracy of derivatives of Galerkin truncated semiflow 11 

3. A-stable Runge-Kutta methods under Galerkin truncation 14 

3.1. Regularity of Galerkin truncated time-discretization 15 

3.2. Convergence of Galerkin truncated time discretization 17 
}J] ■ 3.3. Accuracy of derivatives of Galerkin truncated time discretization 19 
5^ , Appendix A. Stability of contraction mappings 22 

Acknowledgments 30 

References 30 



1. Introduction 

We study semilinear evolution equations 

dtU ^ F{U) :^ AU + B{U) (1.1) 

posed on a Hilbert space 3^ under spectral spatial Galerkin truncation and temporal 
discretization by a large class of A-stable Runge-Kutta methods. The methods 
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2 M. OLIVER AND C. WULFF 

considered permit a well-defined temporal semi-discretization [9, 12]; particular 
examples are Gauss-Legendre Runge-Kutta schemes. The linear operator A of 
(1.1) is assumed to be normal and to generate a strongly continuous, not necessary 
analytic semigroup; _B is a bounded nonlinear operator on y. (This setting includes, 
without loss of generality, cases where A is normal up to a bounded perturbation 
as a bounded non-normal part can always be included into the operator B). The 
examples we have in mind are semilinear Hamiltonian evolution equations such as 
the semilinear wave equation or the nonlinear Schrodinger equation, though for the 
results in this paper we do not assume a Hamiltonian structure. 

Differentiation of the semiflow in time results in multiplication with the un- 
bounded operator A. Hence, in general, the time derivative of the semiflow is only 
well-defined when considered as a map from a subset of D{A) to y [13]; to be 
able to differentiate repeatedly, we assume that B is C^^*^ as map from some open 
set Vk C yk ^ D{A'') to yk for fc == 0, . . . , i^ and iV > K. We prove that the 
semiflow of the truncated evolution equation and its temporal discretization are of 
class C^ jointly in time (resp. stepsize) and in the initial data when considered as 
a map from Vk C 3^k to y, with uniform bounds in the spatial resolution. Analo- 
gous results hold true for the full semiflow and its time-semidiscretization [12]. We 
prove convergence of the space-time discretization on open sets of initial data with- 
out the need of a Courant condition that couples spatial and temporal resolution. 
We then provide estimates on the truncation error of the Galerkin approximation 
of the semiflow, the temporal discretization and their derivatives, and study the 
dependence of the order of the truncation error on the smoothness of the initial 
data. 

We present our results in two versions: we label results that provide uniformity 
of the time interval of existence (for the semiflow) and the maximum time step (for 
the numerical method) on sufficiently small balls of initial data as "local version." 
Assuming more regularity for the initial data, we also obtain results which are 
uniform on bounded open sets so long as B is well-defined and bounded. We will 
label results of this type by "uniform version." 

The standard requirement for the existence of a semiflow is Lipshitz continuity 
of the nonlinearity B [13]. It holds true for a large class of evolution equations and 
will be referred to as condition (BO) in the main text of the paper. Our assumption 
above, made precise as condition (Bl) in the main text, implies Lipshitz continuity. 
Whether the stronger condition (Bl) holds true with K > depends nontrivially 
on the evolution equation and its boundary conditions. It is satisfied by our main 
examples, the semilinear wave equation and nonlinear Schrodinger equation with 
smooth nonlinearities and periodic or homogeneous Neumann boundary conditions. 
It is also satisfied for homogeneous Dirichlet conditions under additional conditions 
on the nonlinearity, see Section 2 below. If (Bl) is not satisfied for sufficiently large 
K, we cannot ensure that the solution U{t) of (1.1) and its numerical approxi- 
mations have enough temporal smoothness to obtain full order convergence of the 
time-discretization independent of the spatial resolution. 

We recall from [13] that the solution to the full semilinear evolution equation 
is obtained as a fixed point of a contraction map, which we consider on the scale 
of Hilbert spaces yo, . . . , yn- Similarly, the Runge-Kutta temporal discretizations 
are functions of the Runge-Kutta stage vectors, which in turn are obtained as fixed 
points of contraction maps. Remaining in this setting, we now consider spatial 
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Galcrkin approximation as a perturbation of these contraction maps. To do so, we 
provide an abstract theory for the stabihty of fixed points under perturbation of 
contraction mappings on a scale of Banach spaces, thereby extending the theory of 
contraction maps on scales of Banach spaces from [12, 16, 18]. This theory provides 
us with a unified framework for the time-continuous and the time-semidiscrete case. 
Let us mention some related results. Spatial spectral Galerkin approximation 
(also called Faedo-Galerkin approximation) is frequently used as a theoretical tool 
for the construction of solutions to partial differential equations; see, e.g., [6, 14]. 
Error estimates for smooth solutions of parabolic problems under spectral and more 
general Galerkin approximations (such as finite element methods) can be found, 
e.g., in [7, 15]. In the parabolic case, there has been a lot of interest in the so-called 
nonlinear Galerkin method which has been shown to have a better convergence 
rate than the standard spectral Galerkin method, see [4, 10] and references therein. 
For analytic initial data, an exponential rate of convergence of the Galerkin ap- 
proximation to the semifiow of the Ginzburg-Landau equation has been shown in 

[5]. 

Hyperbolic problems, namely the semilinear wave equation, and their discretiza- 
tions have been studied, e.g., by Baker et al. [2]. They provide estimates for the 
order of convergence of the spatial Galerkin method of the semifiow for smooth 
enough data and globally Lipshitz nonlinearities under an assumption on the ellip- 
tic projection of the solution; they also treat explicit multistep time discretizations 
of the spatial approximation under a Courant condition that couples the accuracy 
of the Galerkin method with the time-stepsize. Bazley [3] shows the convergence 
of the Faedo-Galerkin approximations of the semilinear wave equation for a spe- 
cial class of nonlinearities on the interval of existence of the continuous solution. 
Verver and Sanz-Serna [17] identify general consistency and stability conditions in 
which convergence of spatial semidiscretizations and of their temporal discretiza- 
tions can be proved. They further verify these conditions for a nonlinear parabolic 
PDE and for the cubic nonlinear Schrodinger equation. In this paper we provide 
a general framework in which those conditions hold true with uniform bounds on 
open sets of initial data. Miklavcic [11] studies a class of parabolic and hyper- 
bolic semilinear evolution equations with a linear part that generates a C'^ semi- 
group, and shows pointwise convergence of the spatial Galerkin approximations of 
the semifiow; he considers nonlinearities B which are Lipshitz on the whole of y. 
Karakashian et al. [8] study a class of implicit Runge-Kutta time-discretizations 
(including Gauss-Legendre methods) and spatial Galerkin approximations for the 
cubic nonlinear Schrodinger equations and prove convergence for smooth solutions 
under mesh conditions that couple spatial and temporal resolution. 

In this paper, the emphasis is on estimates for the spatial Galerkin truncation 
error of the joint higher order derivatives in time and in the initial data both of 
the semifiow and of its temporal discretization. In contrast to [2, 8], our estimates 
for the numerical method hold uniformly in the time-stepsize and do not require 
conditions that couple the spatial and temporal accuracy of the discretization. 
Our results include statements on the pointwise convergence of Galerkin spatial 
semi-discretizations for non-smooth solutions of (1.1) on their interval of existence, 
see Theorem 2.3. These are similar to the results of [3, 11], but include more 
general evolution equations. Our results yield algebraic orders of the Galerkin 
truncation error for smooth, but non-analytic initial data. However, using the 
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methods developed here, it is also possible to obtain exponential estimates for 
analytic data as in [5]. 

The paper is organized as follows. In Section 2, we introduce the class of semi- 
linear evolution equations considered, and show how the semilinear wave equation 
and the nonlinear Schrodinger equation fit into this framework for different types 
of boundary conditions. In this setting, we study regularity and stability under 
Galerkin truncation of the semifiow. In Section 3, we apply a class of A-stable 
Runge-Kutta methods to the semifiow of the Galerkin truncated evolution equa- 
tion and prove results on regularity and stability under Galerkin truncation for 
the temporal discretization which are analogous to the semifiow. We also study 
convergence of the space-time discretization. 

In the appendix, we present a number of technical results on stability of fixed 
points of contraction maps on scales of Banach spaces which are needed in the main 
body of the paper. 

2. Semilinear evolution equations under Galerkin truncation 

We begin by introducing the class of semilinear evolution equations which we 
study in this paper. We then prove regularity of the Galerkin truncated semifiow 
with uniform bounds in its spatial resolution and analyze the dependency of the 
truncation error of the semifiow and its derivatives on the smoothness of the initial 
data. 

2.1. General setting. We consider the semilinear evolution equation (1.1) on a 
Hilbert space y and assume the following. 

(A) A is a normal operator on a Hilbert space y which generates a C°-semigroup 

(BO) B:V^yis Lipshitz. 
Recall that an operator A is normal if it is closed and AA* — A* A. For a defini- 
tion of strongly continuous semigroups (C'^-scmigroups), see [13]. Assumption (A) 
implies that there exists w G K. such that 

Re(spec yl) < u; and ||e*^|| < e"* (2.1) 

for all t > 0. In case A = An + At, where A^ satisfies (A) and At, is bounded, we 
can redefine B a,s B + A\^ and and A as An, whence conditions (A) and (BO) hold 
true. This situation is typical for semilinear wave equations, see Example 2.1. 

For fixed T > and C/° e P let T4^ G C([0, 1]; J) satisfy the fixed point equation 
W = Il{W; U°,T) where, for t G [0, 1], 

n{W; C/°, T){t) = e^^-^t/o + T T e^^^")^-^ B{W{a)) da . (2.2) 

Jo 

When T is small enough, 11 is a contraction on the space Cb([0, 1];3^) so that the 

contraction mapping theorem implies the existence of a fixed point [13]. We then 

define the semifiow <I> of (1.1) by $'^^(C/°) = W{U°, T){t). We sometimes write $* 

to denote the map <&( • , t). 

It is apparent from (2.2) with B — that the £-th time derivative of U{t) is in 

y only if U'^ G D{A^). Hence, we work on a hierarchy of Hilbert spaces defined as 

follows. We set yo = y and, for £ G N, wc define y^ = D(A^) endowed with scalar 

product 

{Ui,U2)y,^{PUi,PU2)y + {\AfQUi,\AfQU2)y. 
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Here P = Pi is the spectral projector of A onto the set {A G spec(j4) : |A| < 1}, and 
Q = 1 — P. This definition of the norm ensures that 

\\A\\y^^^^y^<l and \\U\\y^<\\U\\y^^^ (2.3) 

foraU t/e3^£+i. 

Finally, we write B e C^{V,y) for some V C y ii B € C'^{V,y) and if its 

derivatives are, in addition, bounded and extend continuously to the boundary. 

Then we can state a condition under which $ defines a semiflow on the scale of 

spaces yo,...,yK- 

(Bl) There exist K G No, N e N with N > K, and a nested sequence of 3^fc- 

bounded and open sets 2?^ C yk such that B e C^^'^{Vk,yk) for k = 

0,...,K. 

We denote the bounds of the maps B: Vk ^ yk and their derivatives by constants 
Mfc, M^, etc., for k = 0,...,K and set M = Mq, M' = M^, and so forth. In 
addition to the domains Vq, . . . ,'Dii defined in (Bl), we will sometimes refer to 
'Pk+i which may be any 3^if+i-bounded and open subset of Vk and define 

Rk+i= sup \\U\\y (2.4) 

ueVK+i 

We now give two examples of PDEs that satisfy assumptions (A) and (Bl). 

Example 2.1 (Functional setting for the semilinear wave equation). For the semi- 
linear wave equation 

dttu = dxxU - f{u) (2.5) 

on 7 = (0, 1) with periodic boundary conditions u(0) = m(1), we set U = {u,v) and 

yi^ne+i{I;R)x'HeiI;R) 

for ^ G N. Here, H^(/;R) denotes the Sobolev space of square integrable functions 
whose first £ weak derivatives are square-integrable. Then the operators A and B 
are given by 

^={dl o) ' ^=(l-IPo)i, and B{U) ^ (^_^^^^y (2.6) 

where Pq is the spectral projector of A to the eigenvalue 0. Note that we have moved 
VqAU into the nonlinearity B as Pq^ is not normal. Then the group generated by 
A is unitary on any yi and A generates a C'^-group on y^. So, assumption (A) 
is satisfied. If the nonlinearity / of the semilinear wave equation (2.5) is, e.g., a 
polynomial, then (Bl) is satisfied for any K and N as Hi is a topological algebra for 
e > 1/2 [1]. More generally, if / G C^{D, K) for some iV G N and D C K open, then 
(Bl) holds for K < N; see, e.g., [12, Theorem 2.12]. The same holds true in the 
case of homogeneous Neumann boundary conditions. For homogeneous Dirichlet 
conditions we must additionally require that /^^■'•'(O) = for < 2j < iiT — 1; the 
same restriction on the nonlinearity applies when A is a nonconstant coefficient 
operator and K < 4, see [12, Section 2.5]. 

Example 2.2 (Functional setting for the nonlinear Schrodinger equation). For the 
nonlinear Schrodinger equation 

idtu^ -dxxU + duV{u,u) (2.7) 
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with periodic boundary conditions on / = (0, 1), wc set U = u and identify 

A^idl and B{U) = -iduV{u,u) . (2.8) 

The Laplacian is diagonal in the Fourier representation with eigenvalues — fc^ where 
A: G Z. Hence A generates a unitary group on the square integrable functions 
C2 s jO,2{I; C) and, more generally, on every Heil] C) with ^ e No- So the operator 
A is normal, and assumption (A) holds trivially. In the notation of the abstract 
functional setting of Section 2.1, we choose yg — H2i+i{I;C). If the potential 
Viu,u) satisfies V € C'^+'^+^{D,R'^) for some open subset D cR^ = C, then, by 
[12, Theorem 2.12], the nonlinearity B defined in (2.8) satisfies assumption (Bl) 
for K < N and, in particular, (BO). 

2.2. Spectral Galerkin truncation and convergence. We now truncate the 
evolution equation (1.1) to an A-invariant subspace (Galerkin subspace) as follows. 
For m e N let Pm be the sequence of spectral projectors of A onto the set {A S 
spec(j4) : |A| < m}. Then, assumption (A) implies that 

hm P„iU = U 

for aWU ey, and 

\\A¥mU\\y<m\\VmU\\y (2.9) 

for m (z N. Functions U Cz yi arc well approximated by their Galerkin projections 
PmU- Indeed, setting Q^ — id— P„, 

\\Qn.U\\y<m-'\\U\\y^. (2.10) 

We now introduce the restricted evolution equation 

(2.11) 

where Bm = Pm-B. We write (/)^([/°) to denote the semiflow of (2.11) with initial 
value m™(0) = P™C/°. 

The following theorem provides well-posedness for the projected system on the 
same interval of time on which a solution to the full equation exists, and convergence 
of solutions. 

Theorem 2.3 (Convergence of the projected system). Under assumptions (A) and 
(BO), let U G C{[0,T],T>) be a mild solution to the semilinear evolution equation 
(1.1) with initial value U{0) = U^ . Then there is m^, G N such that for every 
m > TO, there exists a solution Um G C([0,T],2?) to the projected system (2.11) 
with initial value Wm(0) = PmU^ ■ Moreover, 

sup \\U{t)^u^{t)\\y^Q (2.12) 

te[o,T] 

as TO, —!> oo. 

Proof. Local existence of a solution Um{t) of (2.11) is obvious since Am is bounded. 
However, we need to show that the interval of existence is at least [0,T]. We note 
that the solution can only cease to exist if Um leaves the domain 2?, so we proceed 
to prove (2.12) directly. Clearly, 

U{t) - um{t) - e*^Q,„C/" + / e(*-^)^ {B{U{s)) - B,n{um{s))) ds . (2.13) 
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Taking the 3^-norm and noting that, by (2.1), there is c > such that ||e*'^||£(j;) < c 
for t e [0,r], we find that 



\\U{t)-u„M\y<c\\QrnU''\\y + c j \\B{U {s)) - B^{u,r,{s))\\y As 

<c\\Q,nU''\\y + cT sup 



se[o,T] 



<c|lQ„f/"|L + cr sup I 
se[o,T] 



{U{s))\\ 


y + 


Jo 


\\B{U{s))- 


BiUrnisMy 


, ds . 


{U{s))\\ 


y + 


cM'^ 


f\\U{s)- 
Jo 


U,n{s)\\yds 


(2.14) 


tions fra{s) 


^n^B{U{s))\\ 


y converges 


point- 



\frnisi) - frnis2)\ < |!Q™(B((7(si)) - ^([/(.Sa))) Hj; < ||B((7(si)) - ^([/(sa)) || j; , 

the sequence is uniformly equicontinuous. Hence, by the Arzela-AscoH theorem, 
fm converges to zero uniformly as ?Ti — > oo. Thus, applying the Gronwall inequality 
to (2.14), we see for any e > there exists a possibly larger to* such that for 
TO > TO*, \\U{t) — Um{t)\\y < £ so loug as Um{t) docs not leave V. Choosing 
e < dist({C/(s) : s G [0, T]}, dV), we conclude that t in this estimate may be chosen 
as large as T. D 

Let V iz y he open. Wc denote by V^^ some open subset of y which satisfies 

V-^ ^{U eV: d[sty{U,dV) > 5}. (2.15) 

Given (5 > and a hierarchy of open sets Vi C yi ioi £ = 0, . . . , L ioi L E N with 
Vq = V, we define V^^ = V-^ as in (2.15) and, for ^ == 1, . . . , i, 

V^^ = {U eVr- disty{U,dVe) > 5} . (2.16) 

Then, by construction, V^^ C V^\ ior i = 1, . . . , L. Moreover, Bf'{U) C Vi for 
all U G V^ and i = 0, . . . , i where, for any Banach space X and X^ G X, we 
write 

B^{X') = {XeX:\\X~X"\\^<R} 

to denote the closed ball of radius R around X". 

The following theorem provides higher order bounds for the Galcrkin approxi- 
mation error of the scmiflow. 

Corollary 2.4 (Convergence of the projected system - higher order error bounds). 
Assume (A) and (Bl). Let S > be such thatV~^_^^ is nonempty. Then there exists 
TO* such that for all to, > m* and every mild solution U G C([0,T]; 2?]^^-^) of the 
semilinear evolution equation (1.1) there exists a solutionum G C([0, T], P^Hj^x+i) 
to the projected system (2.11) with initial value Um{0) = Pm^° such that 

sup \\U(t)-u^{t)\\y^O{m-''-') (2.17) 

te[a,T] 

The order constants in (2.17) depend only on the bounds afforded by (Bl), (2.1), 
and (2.4), on d, and on T. 
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Proof. As in the proof of Theorem 2.3, we begin with (2.13). Here, we apply Pm 
and rearrange terms to obtain the estimate 

\\U{t)-Urnmy< \\QmU{t)\\y + C f \\B{U (s)) - B{um{s))\\y ds . (2.18) 

Jo 

Due to (2.10), IIQmC^lOllj' ^ Rk+iitt-^^^^- The mean value theorem applies so 
long as Um{s) G V. Then, by the Gronwall lemma as before, we find that (2.17) 
holds true for all m > m^, where choose tti* such that \\U{t) — Um{t)\\y < S for 
t G [0, T] and m > m^ so that indeed Um{s) e T) for s G [0, T] and m > m^. D 

Note that Corollary 2.4 with y replaced by any of the 3^i, . . . , yx readily implies 
that supjg[o,T] 11^(0 ~ Um{t)\\y. = 0{m~^~^~^). However, as B is not assumed to 
map from an open subset of yx+i to ^k+i. Theorem 2.3 as stated does not apply 
with y replaced by yx+i- However, we can still prove the following. 

Corollary 2.5. Under the assumptions of Corollary 2.4, the following is true. 

(a) If N > K + I, we have sup^^rQ j^i ||t/(t) — Um(i)||3^/f+i —5- as to —> oo. 

(b) If N — K + 1, there exists C > such that sup^^^iQ rpA\um(t)\\yj^^-^ < C . The 
bound C depends only on the bounds afforded by (Bl), (2.1), and (2.4), on 
6, and on T . 

Proof. We may assume without loss of generality, that K = {]. (Otherwise replace 
y with yx) Suppose first that N > K + I. Then Theorem 2.3 applies to the 
system of evolution equations 

U^AU + B{U), W^AW + B'{U)W 

with initial value 1^(0) = AU^ + S(C/0) so that W{t) = U'{t). Hence, 

sup \\U'{t)-u'^(t)\\y^O 
te[o.T] 

as m ^^ CO. Since sup^grg yi||i3([/(i)) — Pm-B(um(0)llv = 0{m^^) by Corollary 2.4 
and AU{t) = U'{t) - B{U{t)), we obtain statement (a). 

To prove statement (b), integrate Wm = Awm + B' {ujn)Wm with initial value 
Wm(0) = ^PmC/'^+Pmi3(um(0)) and apply a standard Gronwall argument as before, 
noting that the 3^-norm oiUm{t) is bounded uniformly in ?7i > to^ by Corollary 2.4. 
Thus, supjgFQ m \\u'^{t)\\y < c for some c > depending only on the bounds afforded 
by (Bl), (2.1), and (2.4), on S, and on T. This, together with the bound of B on 
T>, proves (b). D 

2.3. Regularity of Galerkin truncated semiflow. We first introduce some 
notation. For Banach spaces X and y, and j G No, we write £^{y,X) to de- 
note the vector space of j-multilinear bounded mappings from y to X; we set 
£^{X) = £^{X,X). For Banach spaces X, y, and Z, and subsets U <Z X, V d y, 
and W C Z, we write 

FGC^-'"^(i^xV;W) 
to denote a continuous, bounded function F: U x V ^ W whose partial Frechet 
derivatives D^Dyi^(X, F) exist, are bounded, and are such that the maps 

{X,Y,Xi,...,X,)^-D'xBi.F{X,Y){Xi,...,X,) (2.19) 

are continuous from U x V x X^ into £^ {y, Z) for i — 0, . . . , to, and j — 0, . . . , n 
and provided the maps (2.19) extend continuously to the boundary. (The latter 
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is important as we will apply the contraction mapping theorem to maps in such 
classes.) In our setting, V will typically be an interval of time. 

The following theorem provides regularity of the Galerkin truncated semiflow 
with bounds uniform in m under conditions (A) and (Bl) analogous to the regularity 
result for the semiflow $ in [12, Theorem 2.4]. 

Theorem 2.6 (Regularity of the Galerkin truncated semiflow, local version). As- 
sume (A) and (Bl). Choose R> such that V^^ ^ and pick U" e 2?^^. Then 
there is T* = T^{R, U'^) > and m^^iR, f7°) £ N such that for m > tti, there exists 
a semiflow (ply^ of (1.1) of class 

</-»€ fl C^^--''\bI^U'')x[Q,T,];BI'-\U'>)). (2.20a) 

j+k<N 
l<k<K 

(j)m and T* depend only on the bounds afforded by (Bl) and (2.1), on R, and on 
J7°. In particular, 

cj,meC^{Bl'J^{U'') X [0,T,];S^(C/")). (2.20b) 

Proof. The proof is an application of Theorem A. 9 (a) on contraction mappings on 
a scale of Banach spaces. We consider 11 from (2.2) and write the corresponding 
contraction map for the semiflow (pm of the projected system as 

IV^{W-U,h)^^^IV{W-J„,U,h). (2.21) 

We replace N from Theorem A.9 by iV - 1, set /i = T, I = (0, T,), A" = yK with 
U = Uk = B^/2(t/°), w = W, and Zj = C([0, 1]; 3^^) with Wj = C([0, 1]; i3^^ ([/")) 
for j = 0, . . . , K. We now show that the contraction maps !!„ satisfy conditions 
(i) and (ii) of Theorem A.9 for some T^{R, U°) > and m > m^R, U°). We first 
show that Ilm maps each Wq, . . . ,Wk into itself. We estimate, using (Bl) and 
(2.1), that 

\\I1^{W;U,T) U%^ < ||e-^^t/0 - U^^ + ||e-^^(P„t/0 - U')\\y^ 

+ c"^ R/2 + T e"^ Mj . (2.22) 

Choosing T* = T^{R,U^) > sufficiently small, the second line of (2.22) can be 
made less than 3i?/4. Moreover, for a possibly smaller value of T,, there exists 
TO* — m^,{R, U^) such that for aU to > to,, T e [0,r*], and r G [0, 1] the first line 
of (2.22) is less than i?/4. Then, the right hand side of (2.22) is less than R which 
proves that !!„ maps back into Wj. Assumption (Bl) and (A) then imply condition 
(i) of Theorem A.9. To show condition (ii) we estimate, noting that A^ > ii' by 
(Bl), that 

\\Dwn„,{W;U,T)h^cMo.iy,y,))<Te^'^M'^, (2.23) 

so that Ilm is a uniform contraction for all to. > m,, U G U, W E Wj, and 
T e I = (0, T*) for every j = 0,. . . ,K with a possibly smaller value of T,. 

Hence, Hm satisfies conditions (i) and (ii) of Theorem A.9 with bounds and 
contraction constants which are uniform in ?Ti > to,, so that Theorem A.9 (a) 
implies that 0m is of class (2.20a). The simplified special case (2.20b) is a direct 
consequence of Lemma A. 2. D 

Theorem 2.6 does not guarantee that m^ and T, can be chosen uniformly over 
V. The following theorem states that such uniformity can be obtained, however. 
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over domains other than baUs at the expense of stepping up on the scale of Hilbert 
spaces. The situation is analogous to that for the semifiow $; see [12, Theorem 2.6]. 

Theorem 2.7 (Regularity of Galerkin truncated semifiow, uniform version). As- 
sume (A) and (Bl). Choose 6 > small enough such that T>~^,^ ^ 0. Then 
there exists T, — T, {5) > and m^ (S) £ N such that for m > to* the semifiow 
{U,t) f-> 0*„(C/) of (2.11) satisfies (2.20a) with bounds which are uniform for all 
U^ G ^K+i with R = S and such that 

j+k<N 
e<k<K+l 

with bounds which are uniform in m > m*. The bounds on (f)m, "^* dnd T, depend 
oniy on the bounds afforded by (BO) rsp. (Bl), (2.1), and (2.4), and on 6. Moreover, 
(j)m maps into T>k and, when N > K + 1, 

0™ e C^+HVk%, X [0, T,];y) (2.24b) 

with corresponding uniform bounds. 

Proof. We continue to work in the setting introduced in the proof of Theorem 2.6. 
Here, we need to verify that the conditions of Theorem A. 9 are satisfied uniformly 
in U*^ e '^K+i fo'' both 11^ and !!„. First, due to (Bl), each of the Hm is well- 
defined as a map from Wj x U x I into Zj for [/° G ^k+i ^^'^ ^^^ ^^'^ required 
regularity. To show that there is to.* (6) such that !!„ maps Wo , • • • , y^K back into 
itself, we apply (2.22) for every C/° G 2?]~-_|^j^. We bound the first line on the right 
of (2.22) by 

< max(rpe*^C/0||3;. + ||e*'4Q„t/0|| )<e-^i?K+i(T+l/TO), (2.25) 

where Rk+i is defined in (2.4) and j = Q, . . . ,K. Inserting this estimate into 
(2.22), we see that we can choose T* > small enough such that limi ■ ',U,T) maps 
B^'iU°) with R = 6 into itself for aU U° G 2?^+i, C/ G W, to > to*, T G [0,r*], 
and j = 0, . . . ,K. Hence, H^ satisfies the conditions of Theorem A. 9 with bounds 
which are uniform in m > m,, T G (0,T*), and C/° G V]^^^. This shows that 
(2.20a) holds uniformly for t/" G P^"^! and to > to* ((5). 
Next, wc show that 

A0„G fl Ci^''\v],^_^,x[0,n];yk-e) (2.26) 

j+k<N 
l<k<K 

with uniform bounds in to > ?7i*. Consider the linear fixed point equation Wm = 
^m{wm;U,T) with 

+ T f c<-^-^^^''-DBmiw^{<7)){wm{a) + Biwm{a))da . 

(2.27) 

Here, as before, Wm{U, T){t) — (f)'^^{PmU). Integrating the right hand side of (2.2) 
by parts, replacing B with Bm and U''' by Pmt^^ we see that the fixed point Wm 
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of Tim satisfies Wm ~ Aw„i in Zj for j = 0, . . . , K. We consider Ilm with U, 
Zj and T as before, and set Wj = C([0, l\;Br^{0)) with r > large enough that 
n„i(-, U,T) maps Wj into itself for m > to*, U <E U, T <E I. Since 0m is of class 
(2.20a), Lemma A. 6 (a) and Lemma A. 8 (a) imply that n„i satisfies the conditions 
of Theorem A. 9 with N replaced by A^ — 2. Therefore Theorem A. 9 (a) applies and 
proves (2.26). 

Moreover, Bm ° (pm is of class (2.26) with uniform bounds for to. > to* due to 
the chain rule. Lemma A. 6 (a) and the fact that (p^ is of class (2.20a) with uniform 
bounds in t/° and m > to,*. We conclude that dt<j>,n = ^'/'m + ^m ° 0m is also of 
class (2.26) with uniform bounds for m > to.*. 

Finally, as both A0m and dt4'm are of class (2.26), Lemma A. 4 implies (2.24a). 
The simplified special case (2.24b) is a direct consequence of Lemma A. 2. D 

2.4. Accuracy of derivatives of Galerkin truncated semiflow. In this sec- 
tion, we consider how the perturbation of the contraction map H introduced by the 
projection of the evolution equation (1.1) onto the subspace Pm^ propagates into 
derivatives of the resulting semifiow. 

As before, we consider a local and a uniform version of each result; the scales 
we use are defined, separately for the two cases, as follows. In the local version, we 
follow the setting of Theorem 2.6, where we consider initial data 

U eU = B'^_^ (L/°) where X = yx ■ (2.28a) 

he semiflows are considered as maps 

$*,0L:^^f(t^°)^:Vj for j = 0,...,K, (2.28b) 

where m > m*((5, C/°) and i?* = R/2. 

In the uniform version, we follow the setting of Theorem 2.7 where we consider 
initial data 

U eU = P^'^i CX = yK+i ■ (2.29a) 

The semiflows are considered as maps 

$*,0^,: V;, ^y, for i = 0,...,i^ + l, (2.29b) 

for some fixed S > where m > m^{5). 

To handle the complexity of these estimates it is useful to define norms on the 
various objects that contain all combinatorially possible orders of differentiation 
and scale rungs subject to certain relevant side constraints. The need to consider 
such norms arises through the implicit nature of the definition of the semifiow and 
the use of the chain rule. Here, any attempt to estimate a particular derivative 
on a particular rung of the scale will produce terms of all intermediate orders of 
differentiation and scale rungs. We therefore estimate all derivatives at once. 

We have to deal with two different types of objects: contraction maps which 
are functions of three arguments whose corresponding norms are denoted ||| • ||| and 
semiflows which are functions of two arguments whose corresponding norms are 
denoted || • ||. 

In our setting, there are two natural global parameters, namely N, the degree 
of differentiability of the nonlinearity, and K, the number of rungs on our scale as 
defined in condition (Bl). Two more characteristic parameters are needed. First, 
the loss index S which forces the image of the map be estimated at least S rungs 
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down the scale. We will see that a loss of S scale rungs translates into 0{in^^)- 
smallness of the perturbation caused by the projector ¥m- Second, a lowest rung 
index L which forces the estimation of the image of the map to occur at least L 
rungs up from the bottom of the scale. This leads us to define a four parameter 
family of norms for functions 11 = n(w; u, ii) mapping Wk+s x U x I to Zk-e, 

\1^\In,k,l,s = ^_^j^^^_g\\^w^i^tJ^^\\c^(w^+s^uxX:£^Zk+s,£HX:Zk-i))) (2-30a) 

L+e<k<K~S 

for 0<L<K — S<N — S. When studying semiflows, we identify w — W, 
u ^ U, ti ^ T,I = (0,r,), Zj = Ci[0,l];y,), U - B^^^^iU^), X ^ yK, and 

Wj = C([0, l]\B^j{ (L/°)), and use H defined by (2.2). We abbreviate 

inN.K.L-m\N,K,L,0^ (2-30b) 

ll|n|IU,K = llin||lw,K,o,o- (2-30C) 

Functions w — w{u, fi) mapping lA xX to Zj are equipped with the three parameter 
family of norms 

\Mn,k,l= ^.max^ \\^i9iM\c^(uyci-£nx;Z,.,)) (2-30d) 

L+t<k<K 

for < L < K < N, where we abbreviate 

\\M\n.k^\Mn,k,o- (2.30e) 

With this notation, a function (u, ^) i-> !!(«, ^) which does not depend on w satisfies 

in|iAr,K,L,5 = I|n||^„5j4;_5^i . (2.31) 

The next pair of results concerns the stability of the semifiow and its derivatives 
under spectral truncation. 

Theorem 2.8 (Projection error for the semifiow, local version). Assume (A) and 
(Bl) and R > be arbitrary and U" G 2?^^. Let T, = T^{R,U") and m, = 
m^,{R, J7°) be as in Theorem 2.6. Then, for every < P < K , 

ll*-0™IU-p-i,K-p = O(m-^) (2.32) 

where the norm in (2.32) is defined with respect to the spaces (2.28). The order 
constants depend only on the bounds afforded by (Bl) and (2.1), on [/", and on R. 

Proof. We apply Theorem A. 9 to obtain a bound on ||<1' — 4>m\\N-p-i,K-p in terms 
of |||n — nm|||Ar_i^if^o,p, with Zk, X etc. specified above. We already verified con- 
ditions (i) and (ii) of Theorem A. 9 in the proof of Theorem 2.6. Thus, in order to 
prove (2.32), it suffices to show that 

|||n - n,„|||^_i ,^o_p - 0{m-^) . (2.33) 

Writing 

(n-n„0(VK;(7,r)(r) -G,„(c/,T)(T) + /™(VK,r)(r) (2.34) 
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with 

G„(C/,T)(T)=Q™e^"^(7, 

we apply Lemma A. 6 to both tcrm.s on the right of (2.34) in different ways. For the 
first term, since G^ docs not depend on W ^ we can apply Lemma A. 6 with H = id, 
S = P„, and v{U,T){t) = w{U,T){t) = e^^^C/, so that there exists ci > such 
that 

ll|Gm|||/i/_i./f.o,_P ~ W^itt-Wn-I-P.K-P — '-'1 IW"i llw-l,/f,0,P ^ G(m j , 

where the first equality is due to (2.31). Here, and further below, we implicitly make 
use of estimate (2.1) on the bound of the linear semigroup and estimate (2.10) on 
the Galerkin remainder. 

To apply Lemma A. 6 to the second term on the right of (2.34), we identify N 
and K there with A^ — 1 and an arbitrary k E P, . . . ,K here. Then, setting H — id, 
S = P„, u = VK e Wk = C([0, 1]; ^^''(t/O)), and 

v{u, T){t) = w{u, T){t) = r re("-")^^B(VK(a)) da , 

Jo 

Lemma A. 6 asserts that there exists C2 > such that 

\\Im\\N-l-P,n-P < C2 |||Qm ||lAr_i,«,^o,P = Oim^'") . 

Then, by Lemma A. 5, with N replaced by A — 1 and S — P, 

III/ III — 0(m^^] 

The constants ci and C2 depend only on the bounds on B from (Bl) and the bounds 
from (2.1). Altogether, this verifies (2.33), thus concludes the proof. D 

Theorem 2.9 (Projection error for the semiflow, uniform version). Assume (A) 
and (Bl) with N > K + 1, and let S > such that T)'^,^ is nonempty. Let T* = 
r*((5) > and ?7i* — m^^S) be as in Theorem, 2.7. Then, for every < P < K + 1, 

ll<i> - '/'^lU-p-i^K+i-p = 0{m-P) , (2.35) 

where the norm in (2.35) is defined with respect to the spaces (2.29). The order 
constants depend only on S, and on the bounds afforded by (Bl), (2.1), and (2.4). 

Proof. First, we show that 

\\<^~^m\\M-P-i,K-p = Oim-P), (2.36) 

for < P < K on the scale {Zj}j=o,...,a:- To this end, define 11 and Hm, U, y^j 
etc., with R — S as in the proof of Theorem 2.8. 

We have already shown in the proof of Theorem 2.7 that conditions (i) and (ii) 
of Theorem A. 9 hold uniformly in ?7i > m^, and U'^ G 2?^ , -^. Moreover, (2.33) 
holds true uniformly in U^ G 2?]~- , -^. This is easily verified by checking that each of 
the estimates in the proof of Theorem 2.8 holds uniformly under the conditions of 
Theorem 2.9. Hence, Theorem A. 9 (b) implies (2.36). 

Next, we apply Theorem A. 9 to obtain a bound on |||H — Tlm\lN-2.K.o,Pj where 
Hm is from (2.27) and H is defined correspondingly. We have shown in the proof 
of Theorem 2.7 that Hm (and hence also H) satisfy the conditions of Theorem A. 9 
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uniformly for m > m^. Estimating each term of the corresponding anafogue to 
(2.34) via Lemma A. 6 and Lemma A. 8, we find that |||n— nm|||w-2./i'.o,p — 0{m^^). 
Then, Theorem A. 9 (b) implies that ||A$ — A(f)m\\N-p-2.K-p = 0{m^^) so that, 
for < P < if , 

II* - 0mlU_p_i,K-p+i.i = 0(m-^) . (2.37) 

Finally, we prove that 

||5t$ - dt(b^m\\N_p_2,K-p = 0(m-^) . (2.38) 

We note that dt{^ — <t>m) = A{^ — (prn) + B o^ — Bm° <Pm, where the required bound 
on the first term on the right has already been established. For the second term, we 
use Lemma A.6 with H = B, S = B^, ^ = t e I ^ (0,T*), u = U° £ U ^ T^k'+i^ 
wiU'^^t) =. $*(t/0), i;(C/0,t) - C(C^°), X - yK+i, Z, - y,, and W^- - V, for 
j = 0, . . . , if . Hence, there exists a constant c\ such that 

||B O $ - B„ O (\)m\\^_2-P,K-P ^ Cl ll|Qm-B|||^_2,K,0.P + Cl II* - '^mlU-2-P,K-P ' 

To estimate |||Qmi?|||Ar-2,ic,o,Pj we apply Lemma A.6 for each k G P,...,K with 
n = id, E = P„, u = W, v{u, h) = w{u, h) = B{W), U replaced by U^ = >V„, N 
replaced by TV — 1, and K replaced by k. Hence, there is some constant C2 such 
that 

IIQmS(W^)llAr-2-P,.-P < ^2 IIIQm IIU^a.K.O.P = 0(m-^) . 

Then, Lemma A. 5 with N replaced by A^ — 1 and S = P implies 

IIIQmi?|IU_2,K,o,p = 0{m-P) . 
Altogether, this proves (2.38). Finally, (2.35) follows from Lemma A. 4 with L = 
due to (2.36), (2.37), and (2.38). D 

3. A-STABLE RUNGE-KUTTA METHODS UNDER GALERKIN TRUNCATION 

In this section, we study a class of A-stable Runge-Kutta methods which are 
well-defined when applied to the semilinear PDE (1.1) under assumptions (A) and 
(Bl). We prove regularity of spectral Galerkin approximations of such methods 
uniformly in the spatial resolution and derive estimates for the approximation error. 
The class of methods we consider is the same as in [12]. 

Applying an s-stage Runge-Kutta method to the semilinear evolution equation 
(1.1), we obtain 

W = U^l + hs^ {AW + B{W)) , (3.1a) 

U^ = U^ + hh^ {AW + B{W)). (3.1b) 

/VKi\ /B(VKi) 

(^y\ w =\ ■■ , B{w) = ■■ 

ywj \b{w')j 

where W^ , . . . , W^ are the stages of the Runge-Kutta method, 

s s 

{^Wy = Y^ a,j W^ , h^W = Y^ bj W^ , 
and A acts diagonally on the stages, i.e., (v4VK)* = AW^ for 2 = 1,. 
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A more suitable form, required later, is achieved by rewriting (3.1a) as 

W = li{W] U, h) = (id -h^^A)-^ {lU + /laS(VK)) (3.2) 

and 

^^{U) ^ S{hA)U + hh^ (id -h^A)-^ B{W{U, h)) , (3.3) 

where S is the so-called stability function 

S{z)^l + zb^ {id-za)-^l. (3.4) 

We now make a number of assumptions on the method and its interaction with 
the linear operator A. First, we assume that the method is A-stable in the sense of 
[9]. Setting C~ — {z eC: Rez<0}, the conditions are as follows. 
(RKl) The stability function (3.4) is bounded with |S(z)| < 1 for all z e C". 
(RK2) The s x s matrices id— za are invertible for all z £ C~. 
Sometimes, we will also assume that a is invertible. Gauss-Legendre Runge-Kutta 
methods satisfy conditions (RKl) and (RK2) with a invertible [12, Lemma 3.6]. 

We now summarize the analytic properties of the operators appearing in (3.2) 
and (3.3), where we use the convention ||M^||jj» = max^^j^UVK^H^;^. Proofs can be 
found in [12, Section 3.2]. 

Lemma 3.1. Assume (RKl), (RK2), and (A). Then there exist h^, > 0, A > 1, 
(T > 0, and cs > 1 such that 

\\{id-hBA)-'\\y^^y^<A, (3.5a) 

||S(/iA)||3;,^^. <l + cr/i<cs, (3.5b) 

\\haA{id-haA)-^\\y,^y, <1 + A, (3.5c) 
for all h G [0, h^]. Moreover, for any £, n, G No, 

{W, h) ^ (id -haAy^W IS a map of class C^-'^^D^I x [0, K]]y') , (3.6a) 

{W, h) ^ haAiiA -haA)-^W is a map of class c{,-'^^(3^| x [0, h^];^) , (3.6b) 

{W, h) -^ h{\A -haA)-^W is a map of class 4"'^+^^ (3^1 x [0, h^];^) , (3.6c) 
and 

[U, h) ^ S{hA)U is a map of class c'^'^\yi x [0, K]]y) . (3.6d) 



3.1. Regularity of Galerkin truncated time-discretization. Let Wm{U''\h) 
denote the stage vector, w-i^{U°, /i) for j = 1, . . . , s its components, and ipl^lU^, h) 
denote the numerical time-/i map obtained by applying an s-stage Runge-Kutta 
method to the projected semilinear evolution equation (2.11) with initial value 
Wm(0) = PmC/". Their regularity, with uniform bounds in the spatial resolution to, 
is stated in the following theorems which, again in local and uniform version, provide 
the analogue to what is known for the time-/i map \1/ in the time-semidiscrete case 
[12, Theorems 3.14 and 2.15]. 

Theorem 3.2 (Regularity of Galerkin truncated numerical method, local version). 
Assume that the semilinear evolution equation (1.1) satisfies conditions (A) and 
(Bl). Apply a Runge-Kutta method'^ subject to conditions (RKl) and (RK2) to it. 
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Choose R>0 such that 2?^^ ^ and pick U° G 2?^^. Let R^ = R/{2 max{cs, A}) 
with A and cs from (3.5). Then there exists to* = m^(R, U'^) and /i, = h^{R, U'^) > 
0, such that for to, > ?Ti, there exists a stage vector Wm whose components w^^ as 
well as the numerical time-h map il^l^{U,h) = tp!^{U) are of class (2.20a) with T, 
replaced by /i,. The bounds on Wm, ipm cind ft,* are independent of m and depend 
only on the bounds afforded by (Bl) and (3.5), on the coefficients of the method, 
R, and U° . 

Proof. As in the proof of Theorem 2.6 on the regularity of the scmiflow (/>„, we 
apply Theorem A. 9 (a) on contraction mappings on a scale of Banach spaces. Here 

we set w = W, Zj = y^, and Wj = B^J (lt/°) for j = 0, . . . , K. We further identify 
H = h,I={Q, K), and X = yK with U =Uk ^ Br,{U") C yx- 

The map 11 for the stage vector W is defined by (3.2); we write the corresponding 
map Iim{W,U,h) = F^Il(W,¥mU,h), analogous to (2.21) for the semiflow. We 
now show that the differentiability assumption on B is such that !!„ satisfies the 
conditions of Theorem A. 9 for to, > vn, with a suitable choice of ?n,. 

First, we show that !!„ maps each Wq, . . . , Wx into itself uniformly for h G 
(0, ft,*) and U <eU. By Lemma 3.1 we estimate, for W £ Wj, 

\\n^{W;U,h) - lU^'Wy. < ||(id-(id-ftaA)-i)l[/0||y. + \\iid-haA)-'lQmU''\\ys 

+ Ai?* + ftA||a||Mj. (3.7) 

Using Lemma 3.1, we can find h*{R,U^) and m*{R,U'^) such that for m > m* 
and ft e (0, ft*) the first line on the right of (3.7) is less than R/A. By possibly 
shrinking ft* further, the second line is less than 3i?/4, so that n„i( ■ ■,U,h) maps 
Wj into itself. Second, assumptions (Bl) and Lemma 3.1 ensure that Ilm satisfies 
condition (i) of Theorem A. 9 for to, > to,*. The contraction estimate, condition (ii) 
of Theorem A. 9, follows from 

\\DwnUW;U,h)\\y^^y^ <h Wild -ha A)-^a\\y^^y^ \\DB^{W)\\y^^y^ 

<ftA||a||Mj (3.8) 

for j ~ 0, . . . , K. Thus, by possibly shrinking ft* again, the right hand bound can 
be made less than 1, and condition (ii) is met for m > to,*. Thus Theorem A. 9 (a) 
implies that w-j^ is of class (2.20a) for j = 1, . . . ,s. The same holds true for ■0™ 
due to the chain rule on scales of Banach spaces, Lemma A. 6 (a), applied to (3.3) 
using Lemma 3.1. D 

As for the semifiow 0m, there is also a uniform version of this result. 

Theorem 3.3 (Regularity of Galerkin truncated time-discretization, uniform ver- 
sion). Assume (A) and (Bl), as well as (RKl) and (RK2). Pick S > such that 
T>]/,i is nonempty. Then there is ft, = ft* ((5) > and to,* ~ m*((5) such that for 
TO, > ?Ti* the statements of Theorem 3.2 hold true with R = S and bounds which are 
uniform in [/" G T^k+i o-nd to > m*. Moreover, for to, > ?Ti*, the components w-j^ 
of the stage vector Wm{U, ft) are of class (2.24a) with T* replaced by ft*, and, if the 
Runge-Kutta matrix a is invertible, the numerical time-h map ipm is also of class 
(2.24a) with T* replaced by ft*. The bounds on Wm, ipm and /i* are independent of 
m > to* and only depend on the bounds afforded by (Bl), (2.4), and (3.5), on the 
coefficients of the method, and on d. 
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Proof. Let Zj = yj, Wj = B§ (If/"), and U = B^'^ (U^) for j = 0, . . . ,K as in the 
proof of Theorem 3.2, taking R = S and R^, as in Theorem 3.2. First, due to (Bl), 
the map !!„ is well-defined from Wj x U x I into Zj with the required regularity 
properties and with bounds that are uniform in m and t/° € T^k+i- ^'^ show that 
Ilm maps Wj back into Wj for m > m*((5) with a suitable choice of m*, note that, 
forj=0,...,K, 

\\{xA-{xA-h^A)-H)U''\\ys + \\{{A-h^A)-^mmU''\\y. 

<h max \\aA{[d-ssiA)-^lU"\\ys + A ||lQ„C/°|Le 
se[o,h] -^3 -^3 

< (/lA^ ||a|| + A/m) sup \\lJ\y, 

< (/iA2||a|| +A/m)i?K+i, (3.9) 

where Rk+i is defined by (2.4). Inserting this estimate into (3.7), we see that we 
can choose h^,{5) > small enough and m^,{d) big enough such that !!„(•;[/, /i) 
maps each yVj into itself and, due to (3.8), such that !!„ is a contraction on each 
Wj uniformly for U^ , U G Vj^,^, and h G [0, ft.*]. So the conditions of Theorem A. 9 
are satisfied uniformly for m > to,, h G (0, ft.,), and [/" G ^^ 1 1. 
Applying A to n,„(VF, [/, ft) yields 

Aw„i ^ AU^{w„i,F^U,h) = (id-ftaA)"HAP„L/ + fta^(lC/-ftaA)"ip„B(w,„). 

Using Lemma A. 6 (a), the chain rule on a scale of Banach spaces, together with 
the estimates of Lemma 3.1, we find that Aw-l^ is of class (2.26) for j — 1, . . . , s. 
Moreover, on the {K + l)-scale Zj — yj for j = 0, . . . ,K and Zk+i = yK with 
Wj = V'j for j = 0, ... X and W^+i = 2?^, and with U ^ X>^+i, X ^ yx+i, and 
I ~ (0, ft,), the map Urn satisfies conditions (i) and (ii) of Theorem A. 9. (Here we 
have used Lemma 3.1, in particular (3.6c), once again.) Therefore, 9/iW^ is of class 
(2.26). Then Lemma A. 4 implies that w^ is of class (2.24a). If a is invcrtible, we 
can use (3.3), (3.6b) and the chain rule in the form of Lemma A. 6 (a) to show that 
V-" is also of class (2.24a). D 

Remark 3.4. When a is not invertible, an according modification of the proof of 
Theorem 3.3 yields the weaker statement 

^™e n Cp+''(2?^^+, x[0,T,];2?) (3.10) 

j+k<N-l 
k<K 

for m > ?Ti, with bounds that depend only on the bounds afforded by (Bl), (2.4), 
and (3.5), on the coefficients of the method, and on S. An analogous statement 
holds true for * [12, Remark 3.22]. 

3.2. Convergence of Galerkin truncated time discretization. Next, we prove 
a convergence result for the time-semidiscrctization of the projected system. The 
error bounds are uniform in the spatial truncation parameter m. 

Theorem 3.5 (Convergence of time discretization of projected system). Assume 
that the semilinear evolution equation (1.1) satisfies condition (A), and apply a 
Runge-Kutta m,ethod of classical order p subject to conditions (RKl) and (RK2) 
to it. Assume further that (Bl) holds with K > p. Pick 6 > such that 'P'p+i is 
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non-empty, and fix T > 0. Then there exist positive constants h^, m^, a, and C2 
which only depend on the bounds afforded by (Bl), (3.5), on the coefficients of the 
method, and on S, such that for every U'^ satisfying 

{$*(t/O):ie[0,T]}cP;+\, (3.11) 

h € [0, ft.*], and m > m^,, the numerical solution ('*/'„)" (^°) lies in D and satisfies 

\\{i^':nnu")-4'7^'{U'')\\y < C2e^^""ftP (3.12) 

so long as nh < T . 

Proof. Convergence of the time seniidiscretization under the above assumptions 
can be proved by a standard Gronwall argument, see [12, Theorem 3.24] where 
condition (Al) there is always satisfied in the setting here; it is stated as (3.5b) in 
Lemma 3.1. Then [12, Theorem 3.24] asserts that whenever a semiflow <I> satisfies 
(3.11), there exist constants ci and C2 which only depend on the bounds afforded 
by (Bl), (3.5), on the coefficients of the method, and on 'DZ,^^ such that 

||(*'')"([/")-$"''(t/0)||y <C2e"i"''/iP (3.13) 

so long as nh < T . 

Here, we need to apply this result with $ replaced by (jjm and ^ by ■(/;„; and show 
that we obtain uniform bounds in to > to,* and [/" satisfying (3.11). So we have to 
show that if condition (3.11) holds, then there is an analogue of this condition for 
the truncated system which is valid for all to, > to,* In other words, wc have to find 
a 3^p+i-boundcd set Pp+i C Vp n 3^p+i and 5 > such that 

{^l{U''):te[0,T]}cV;l, (3.14) 

holds for all f7° satisfying (3.11) and all to, > m*. Applying Theorem 2.4 with 3^p 
in place of y and 2?p in place of 2?, wc find that there is some m^ G N such that 

sup \\^\U')~<l>UU")\\y^<S/2 



and 



sup UUU')\\y^^,<C 

te[o,T] 



for some C > 0, all to, > tti*, and all J7° satisfying (3.11). Thus, with 6 = 6/2 and 

Pp+i=2?pnintS^;+/(0), 

where int(W) denotes the interior of a set W of a Banach space X, condition (3.14) 
is satisfied for all m > to,*. This completes the proof. D 

By combining this theorem with Theorem 2.4 we obtain convergence of the space- 
time discretization to the semiflow $*(C/°) of order 0(/iP) + 0(m^^^^) fort e [0,r] 
and m > m^ with uniform bounds for all C/" satisfying (3.11). In particular, we do 
not require a coupling between spatial resolution m and temporal resolution h for 
this convergence result. 
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3.3. Accuracy of derivatives of Galerkin truncated time discretization. 

Results corresponding to Theorems 3.2 and 2.8 hold true for the stability under 
spectral truncation of the numerical stage vector and its derivatives. 

Theorem 3.6 (Projection error for the numerical method, local version). Assume 
(A), (Bl), (RKl), and (RK2). Fix R > such that D^ is nonempty and choose 
U° e Vj^^. Let K = h4R,U") > and m, = m,^{R,U^) he as in Theorem 3.2. 
Then for every < P < K , 

W-w^\\j,_^_PK_P = 0(m-P) (3.15a) 



II*'' - ^t\\N-i-p.K-p = 0{m-P) , (3.15b) 

where the norm in (3.15) is defined with respect to the spaces (2.28). The order 
constants depend only on the bounds afforded by (Bl) and (3.5), on the coefficients 
of the method, and on R. 

Proof. The proof of (3.15) is an application of Theorem A. 9 on the stability of 
contraction mappings where, as in the setting of Theorem 3.2, 11 is defined by 
(3.2), Em == Pm o n, and we set 

Wj=bI'{1U^) where Z^^y^ (3.16) 

for j — 0, . . . , K. We further identify I ~ (0, /i,), w ~ W, fi ~ h, and X = yx 
withl{ = llK=Bl'^{U°). 

We already verified in the proof of Theorem 3.2 that conditions (i) and (ii) of 
Theorem A. 9 hold with uniform bounds for m > m^. Thus, Theorem A. 9 (b) yields 
a bound of the form (3.15a) provided we can show that 

|||n - n,„|||^_i ,^0 p = Oim-P) . (3.17) 

Writing 

n(VK; U, h) - n„(VF; U, h) = (id -hs^A)-^ IQmt/ + /la (id -/laA)"! <^^B{W) , 

(3.18) 

we apply Lemma A. 6 to both terms on the right as follows. For the first term, 

we take II = id, E = P^, and v{U, h) — w{U, h) = (id —hzA)^^lU to conclude that 
there exists ci such that 

lll'-'m|lw-l,K,0,P ^ \\^rn\\N-P-l.K-P — '^i\\\^rn\\\N-i^K,O.P ~ ^('^ ) i (3.19) 

where the first equality is due to (2.31), we recall Lemma 3.1 for the differentiability 
properties of (id— /laA)^^, and note that the final statement is due to estimate 
(2.10) on the Galerkin remainder. 

To estimate the second term on the right of (3.18), we apply Lemma A. 6 for 
each K G P,...,K with II = id, S = Pm, u = W, and v{u,h) — w{u,h) ~ 
ha {id-haA)-^B{W), U replaced by U^ = B|~(1C/0), N replaced by TV - 1, and 
K replaced by k. Hence, by Lemma 3.1, there is some C2 such that 

|1Q„ ha {Sd-haA)-^B{W)\\^_^_p^^_p < c^ |||Q™||U_i,,,o,p = 0(m-^) . 
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Then, Lemma A. 5 with N replaced by A^ — 1 and S = P iniphes 

IIQ™ ha (id -haA)-^B\\\j^_^j^,^p = 0(m-^) . (3.20) 

The constants ci and C2 depend only on the bounds from (Bl) and the bounds from 
Lemma 3.1. Altogether, we have proved (3.17); the proof of (3.15a) is complete. 
To prove estimate (3.15b), note that by (3.3), 

*"(C^) - i^t(U) = S{hA)Q„,U + b^ {J{W{U, h),h) - F^JiwUU, h),h) , (3.21) 

where 

J{W] U, h) = h{\A -^haA)-^B{W) , (3.22) 

so that 

+ s\\h\\ \\JoW-¥^Jow^\\j^_^_pj^_p. (3.23) 

We estimate the first term of (3.23) using Lemma A. 6 with 11 = id, E = P„i, u ~ U, 
Zj = yj, W-j = Bp'{U") for j = 0,...,K, and wiu, h) = v{u, h) = S{hA)U . For 
the second term of (3.23), we use Lemma A. 6 with IT = J, S = Pm>/, Zj and 
Wj from (3.16) as before, w{U,h) = W{U,h), and v[U,h) = Wrn{U,h). Thus, by 
Lemma 3.1, there exists C4 such that 

II*'' - ^^t\\N-l-P^K-P < C4 ll|Qmt/|||jv_i,K,o,P + C4 IIIQm J|IU_i^k^o,P 

+ C4,\\W -w„,\\j^_i_P^K-p- (3-24) 

The first term is 0{m^^) by (2.10); to obtain the required estimate for the second 
term wc proceed as in the computation proving (3.20), but with /i(id —haA)^^B in 
place of /ia(id — /iaA)"^i?; the third term is 0(m^^) by (3.15a). D 

Theorem 3.7 (Projection error for the numerical method, uniform version). As- 
sume (A), (Bl), (RKl), and (RK2). Choose S > small enough such that 2?x+i 
is nonempty, and let /i* = h^{6) > and m* = to* (J) he as in Theorem 3.3. Then 
(3.15) holds true with respect to the uniform setting (2.29). Moreover, for every 
Q<P<K + 1 andN > K + 1, 

\\W - u;™lU_i_p,^+i_p = Oim-P) (3.25a) 

and, for a invertible, 

11*'^ - i't\\N-i-p,K+i-p = 0{m~P) , (3.25b) 

where the norm in (3.25) is defined with respect to the spaces (2.29). The order 
constants depend only on the bounds afforded by (Bl), (2.4), and (3.5), on the 
coefficients of the method, and on 6. 

Proof For each [7° e P^-^^^, define H, n„, Wj = sJ(lC/0), Z^ = 3^|, I = (0, h,), 
etc., as in the proof of Theorem 3.2. We first show that (3.15) holds true with respect 
to the uniform setting (2.29). We already verified in the proof of Theorem 3.3 that 
conditions (i) and (ii) of Theorem A. 9 hold uniformly in t/'^ G 'Dj^,^ and m > to*. 
Further, by checking uniformity of all required estimates, we verify that (3.17) in 
the proof of Theorem 3.6 hold uniformly in U^ G 2?^ , j^. Applying Theorem A. 9 
(b) for all U^ G T'k^i then implies that (3.15a) holds in the uniform setting (2.29). 
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Similarly, (3.24) holds uniformly in t/° G '^k+i ^^ that (3.15b) holds with respect 
to the uniform setting (2.29). 
We next show that, for N > K, 

\\AW - Awm\\N_i_p^K-p = 0{m-^) (3.26a) 

and, for a invertible, 

\\A^'^ - A^t\\N-i-p.K-p = 0{m-P) , (3.26b) 

in the uniform setting (2.29). We apply A onto (3.18) as well as onto the expression 
for '^^ — ^m- The resulting difference expressions are then estimated as in the proof 
of Theorem 3.6 using Lemma 3.1, in particular (3.5c). Now, we aim to show that, 
for A^ > X + 1 and < P < X, 

119^14^ - dhWm\\i^_p_2,K-p = 0{m^^) (3.27a) 

and 

119/.* - dhyJm\\N-P-2,K-p = 0{m-'') . (3.27b) 

To prove (3.27a), we apply Theorem A. 9 on the stability of contraction mappings to 
the pair 11 from (3.2) and 11^ = PmoH, but this time on the (K + l) scale Zj — y!- 
iorj = 0,...,K and Zk+i = y^ with W^ = P^ for j = 0, . . . , if and Wk+i = ^k- 
Set, as before, U — X'^^-^, X = 3^k+i, and I = (0, /i,). Due to (3.6c) and (Bl), 
the map 11^ satisfies the assumptions of Theorem A. 9 in this setting for m > to*. 
We obtain |||n — Ilm\\\N-i,K+i.o,p — 0{m~^) as in the proof of Theorem 3.6. By 
Theorem A. 9 (b), this implies ||VF — Wmlljv-p-i k+i-p = 0{m^^) with respect to 
the above defined hierarchy, and in particular (3.27a). 

Estimate (3.27b) is proved similarly. We estimate the norms of both terms in 
(3.21), with J as in (3.22). First, as in the proof of Theorem 3.6, using Lemma A. 6 
with n == id, E ^ P„, u == U, Zj = y^, Wj ^ Br' ([/") ioT j ^0,...,K+ 1, where 



^y 



r > is such that Vk+i C £?r'^^^(0), and w{u,h) — v{u,h) = S{hA)U , we obtain 

\\S{hAn^U\\^_P_^j,^^_P = 0{m-P). 
Thus, in particular, 

\\dHS{hA)Q„,U\\j,_p_2^j,_P - 0(to-^) (3.28) 

holds with respect to the uniform setting (2.29b). 

We now estimate the second term of (3.21). Consider the [K + l)-scale from 
above, i.e., Z, = y^ for .?■ = 0, . . . , if , Z^+i = y^. W, = 2?| for j = 0, . . . , if, 

Wk+i = Vk, U = ^xti' -^ = ^^+i' ^"d I = (0,/i,). Due to (3.6c) and (Bl), 
the maps J and Pm J satisfy the assumptions of Lemma A. 6 in this setting, and, as 
we have seen above, W and Wm also satisfy the conditions of Lemma A. 6 for this 
choice of scale and m > to*. We obtain \\\QmJ\lN~i.K+i.o,p = 0{ra^^) as in the 
proof of Theorem 3.6, and 

\\J oW - ¥raJ O Wra\\M^l-P^K+l~P = 0{m^^) 

for the above choice of scale. This implies that, with respect to the uniform setting 
(2.29b), we have 

\\dh{J oW - P„iJ o Wm)\\N-2-p,K-p = 0{m^^) 
which, together with (3.28) and (3.21), implies (3.27b). 
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Finally, Lemma A. 4, given that (3.15) holds in the uniform setting (2.29b) and 
together with estimates (3.26) and (3.27), implies (3.25). D 

Remark 3.8. If, in the setting of Theorem 3.7, the matrix a is not assumed to be 
invertible, then (3.27b) still holds, of. Remark 3.4. 

Appendix A. Stability of contraction mappings 

Abstract contraction mapping theorems on a scale of Banach spaces have been 
obtained in [12, 16, 18]. For the results in this paper, we must, in addition, estimate 
the stability of the fixed point under perturbation of the contraction map. 

For K E Nq, \ct Z = Zq D Zi D . . . D Zk be a scale of Banach spaces, each 
continuously embedded in its predecessor, and let Vj,Wj C Zj be nested sequences 
of sets. Let Af be a Banach space, and let U C X and I C M be open. We note 
that all results in this section easily extend to the case where I is an open subset 
of RP. We may assume that HwH^r^ < ||w||2j+i for aU w € Zj+i. (If this is not the 
case, we inductively equip Zj^i with the equivalent norm || • Wzj+i + \\ ■ \\zj-) 

As detailed in Section 2.3, we use the following additional integer indices. The 
minimal regularity we guarantee for the image space of the function considered is 
scale rung L, the "loss index" S indicates how many rungs on the scale the range 
of a function is down relative to its domain, and N denotes the maximal regularity 
of the function. We assume Q<L<K — S<N — S. We work with the family of 
spaces 

Civ,K,L,s({Vj},Z^,Z;{>Vj})= fl ci-'^''^(Vfc+sxiYxZ;Wfe_,), (A.la) 

i+j+k<N-S 
L+i<k<K-S 

endowed with norm (2.30a), and abbreviate 

Civ,/f,L({V,},W,Z;{>V,}) =CAr,x,L,o({V,},W,I;{W,}), (A.lb) 

Cjv,k({V,},W,Z;{W,}) -CAr,KAo({V,},i^,I;{>V,}) (A.lc) 

with corresponding norms (2.30b) and (2.30c), respectively. We note that any 
function of class (A.la) has a maximal number oi N — L — S derivatives in its first 
and second argument on the lowest admissible domain scale Z^-^-s- 
Furthermore, let 

Cw,k,l(^,I;{Wj})= n Cl^-'\U X I;Wk-i) , (A.ld) 

j+k<N 
L+l<k<K 

endowed with norm (2.30d), and abbreviate 

Cn.k{U,I; {W,}) = Cw,K,o(W,2:; {W,}) (A.lc) 

with corresponding norm (2.30e). For future reference, we note the following. 

Remark A.l. When a map 11 G Cn.k.l.s({Vj},U,I;{Wj}) does not depend on 
w, it can be interpreted as an element from Cis[,k.l{U,I] {VVj}) where, by (2.31), 

|||n|||7V,K,L,S = \\^\\n-S,K-S,L- 

We simply write Cn.k,l,s and Cm.k.l when the arguments are unambiguous. We 
also write 

d^JV{w{u,^i);u,^i) ^ d^Jl{w;u,ii)\^^_^^_^^^^ = (a^B o w;)(m,^) 



GALERKIN APPROXIMATIONS OF A-STABLE RUNGE-KUTTA DISCRETIZATIONS 23 

to denote partial /i-derivatives vs. D^(n(w(u, /i), u, /i)) to denote full /i-derivatives. 
We begin with four technical lemmas which can be proved by simple index arith- 
metic. Details can be found in [12]. 

Lemma A. 2. If N > K then, with W = Wq, 

Cw^KiU,!; {W,}) C Cf (W X Z; W) . 

Lemma A. 3. Suppose that w e Cn,k,l{1^,I', {VVj}) and that the map {u,u,^) h^ 
Duw{u, ii)u is of classCM.K.L{UxBi{Q)^I]{Zj]). Thenw G CAr+i,_fs-,L(^/,I; {VVj}) 
and 

W^Wn+I.K.L^ sup |iD„WM||jY^^+ ||w7||jY^^. 

■ ■ ll"IU<i 

Lemma A.4. When N > K, w e Civ,K+i,L+i(W,I; {Wj}) n Cjv,l,l(W,2:; {W^}), 
and d^w e Cn-i,k,l(1^,I] {Zj}), then 

weCN,K+i,L{l^,X-{Wj}) 

and 

M\n,K+IX < \\M\n^K+IX + 1 + \\M\nXX + \\dM\N-l,K,L ■ 

Lemma A. 5. We have 

S<K<K 



IPIIIc„.K.L,s({V,},W;I;{W,}) '^ 5™^<v"^"c„-s,»-s,i(V„ X W,I;{W, }) ' 

where ~ denotes that left hand and right hand sides provide equivalent norms on 
Cn,k,l,s- 

We now prove a stability result for fixed points of contraction mappings, i.e., we 
want to bound norms of differences of fixed points in terms of norms of differences 
of contraction maps. To do so, we first need to look at a corresponding stability 
result for compositions. 

The following lemma states that the difference between two functions which are 
both compositions of functions can be estimated by the difference of the outer 
functions and the difference of the inner functions, and that the same holds for 
derivatives of the difference. Here S is the minimal smoothness of the image of the 
inner functions and of the domain of the outer functions, K is the number of scales 
and A^ — S" is the maximal number of derivatives of the inner functions. Finally, P 
is used to relax the required smoothness in the estimates. 

Lemma A. 6 (Stability of compositions). For 0<S + P<K<N, suppose 
If = n(w; u, fl), S — S(w; u, fx), w — w{u, fx), and v — v{u, /i) satisfy 

(i) n,^ eCN+i.K.o.s{{y^j}M,i:{z,}); 
(ii) w,t;GCjv,if,s(Z^,2:;{Wj}). 

Then: 

(a) flow G CN-S,K-siM,I]{y^j}) o-nd |lnotf;||^_5 ;^_5^^ can be bounded by a 
polynomial with non-negative coefficients in HHWn.k.l.s (ind \\w\\n^k,s+l; 
the same holds true for Sow. 
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(b) There is some c > which is a polynomial with non-negative coefficients in 

|||n|||jv+i,A:,o,s, ||w||Ar,A:,S; |||S|||Ar+i,_R-,o,s, and \\v\\n,k,s such that 

\\U0W-^0v\\j^_p_gj^_p_g < C\lU-m^j^^„p^g+c\\W-V\\^_pj^_pg. (A.2) 

Remark A. 7. Part (a) was already shown in [12, Lemma A. 6] and is the chain rule 
on the scale of Banach spaces. It will be our main tool for obtaining estimates on 
the scale of Banach spaces for compositions of maps of the form (11 o w) {u, /x) = 
Il(w{u, fj,); u, fi). The essence of the result is very natural: When the outer function 
n loses S rungs on the scale, the inner function w must have minimal regularity 
L = S and the composition maps at best into scale rung K — S. 

Proof. We prove part (b) only. It follows the same pattern as the proof of part (a) 
with the additional difficulty that we need to carefully keep track of differences in 
the various spaces. We proceed by induction in N and K. For N ~ K = P + S, 

\\U0W~^0 V\\^^^^^,^^ < ||n O u; - n O l'||c(;^xI;Z) + P O i; - S O w||c(;^xI;Z) 
< Co ||ui - ^'|lc(iyxI;Zs) '^ 111"'^ ^ ^ip+S,P+S',0,P+S' ' 

where, by the mean value theorem, cq = |||n|||p+s+i,p+s,o,p+s- 

Let us now increment N holding P and K fixed. Let B = i3j^(0). By Lemma A. 3, 
it is sufficient to derive the claimed upper bound for the Cn-p-s.k-p-s{^ x 
B,I; {yVj}) norm of the function which maps {{u, w), /i) E {U x B) x I to 

D„(n o w — T, o v) u = (duTlu) ow— (9„Sm) o v 

+ tl o (Duw n) - Y. o (DuV n) , (A.3) 

where II and E are defined in (A. 8) below. 

To estimate the first line of the right of (A.3), we define Ili{w;{u,u),ii) = 
duTl{w;u, fi) u and Si(w; (■u,u),/z) = du^(w;u, fi) u. Then, by induction hypothe- 
sis, there is a constant ci which is a polynomial in |||n|||jv+i,x.o,s ^ |||ni|||Ar^i<-_o.Sj 
\m\N+hKfi,s, \\w\\n+i,k,s > \\M\n,k,s^ and ||w|U+i,k,s such that 

llHi o w; - El o v\\^_p_sj^_p_s 

< Ci lllni - ^l\lN^Kfi,P+S + Cl Ik - v\\n-P,K-P,S 

— Cl iin — ^\\\N+l,Kfi,P+S + "^1 11^ ^ ^\\n-P,K-P,S ■ 

To estimate the second line of the right of (A.3), fix 

r = max{||w||jY^i^^g, M^^^j^g} (A.4) 

and set Vj — Br^ {0) for j ~ S, . . . ,K. By Lemma A. 8 (a) below, the maps 11, E 
satisfy condition (i), i.e., 

n,±eCN+uKfi.sm},U,I;{Zj}). 

For fixed u Q B, a direct estimate verifies that w((u,w),/i) = T>^w{u,pl)u and 
'0((w, u), /J.) = D„w(u, /i) u map {U x B) x I into each of the domains Vs, ■ ■ ■ , Vk of 
n and E. Hence, w and v satisfy assumption (ii). Then, by induction hypothesis. 
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there is some constant C2 such that 

lln o u; - S o w\\n-P-S,K-P-S 

< C2 |||n - S|||w,K,O.P+S + C2 \\w - v\\n-p,k-p,s 

< C2 |||n - ^lN^Kfi,P+S + C2 \\w - v\\n+1-P,K-P,S ■ (^-5) 

Wc note that C2 is a polynomial in ||io||Ar,K,5, ||w||Af,K,5, |||n|||7v+i,if,o,s which, by 
Lemma A. 8 (a), is bounded by a polynomial in ||w||Ar+i,K,Si |||n|||7v+2,A:,o,S7 and r, 
and |||S|||7v+i,if,o.s which is bounded by a polynomial in ||w||Ar+i.K.Si 111^111^+2.^.0,57 
and r. 

To estimate the term |||n — 'S.\In.k.o.p+s in the last inequality above, note that 
the maps 

n2(w; (w, u), /i) = D^n(iy; u, /i) w and S2(w; (w, u), /i) = Dt„S(w; u, /i) w 

satisfy, ioi P + S < k< K, 

n2,S2 eCAr+i,«,o,s({Wj},V. xZ^,Z;{Z,}). 

So the induction hypothesis applies once again, asserting that there is a constant 
C3 such that 

llHa o w - E2 o v\\j^_p_g^^_p_g < C3 |||n2 - ^2\lN,^fl,P+S + C3 \\w - v\\j^_p^^_pg 

where, for all P + 5* < k < K, 

1112 - ^2\In^i^,o,P+S - '"' in ^ '^tN+l,Kfi,P+S 

and C3 is a polynomial in |||n|||Ar+2,K,o,s, |||S|||Ar+2,K,o.5, \\w\\n+i,k, \\v\\n+i.k, and 
r. Therefore, by Lemma A. 5, there is some constant C4, which is independent of 
112, S2, f, and w, such that 

ll|n - A\N,Kfi,P+S ^ Il|n2 o w - E2 o w|||m^^^o.P+5 

- ^^ P+™<?<K '"^^ ° W' - ^2 o f'||lAr_p_5,«_p_s 

< r C3 C4 |||n - S|||Ar+i,_R-,0,P+5 + C3 C4 ||w - « IU-P,K-P,S ' 

This last assertion is summarized in Lemma A. 8 (with N replaced by A^ — 1). This 
concludes the inductive step in N . 

Second, we prove that the conclusion also holds when we increment K — S when 
K < N, holding N fixed. By Lemma A. 4, 

||n o u; - S o I'lU^p^s.K-P-S+l < ||n o u; - S o "ylU-P-S.K-P-S+l,! 

+ ||n o w - S o w||^_p_5^o + llD^Cn o W - S o w)||Ar_p_5_i,^_p_5 • (A.6) 

To estimate the first term on the right, note that we can apply the induction 
hypothesis on the translated scale Zj = Zj+i, Wj = Wj+i- Thus, there is a 
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constant C5 with the required polynomial dependence such that 

iinoui — SowiL ii/T!-^\\ 

" iiCjv-p-s,K-p-s+i,i(W,I;{2j}) 

= llnow- Sow|| ajT S-7 l^ 

= C5 |||n - S|||g^^^^^_^^^(^y^^j^^ 2-._j2^|) + C5 llw - f'llc„_p,K+i-p,i + s(W,X:{23}) ■ 

For the second term on the right of (A. 6), we apply the induction hypothesis on 
the trivial scale, obtaining that there is a constant cg with the required polynomial 
dependence such that 

||n o w - S o -ylljv-p-s^o ^ C6 |||n - S|||jY_p_^5_o p_^5 + C6 ||u - u^lU.p^s^s ■ 
For the third term on the right of (A. 6), we estimate 

+ ||n O d^W -to df,v\\^_p_g_^j^_p_g . (A.7) 

To estimate the first term on the right of (A.7), notice that H, S G Cn+i.k+i.o,s 
implies 9^11, 9^S S Cjv+i,_r:+i,o,5+i- Since w,v £ Cn,k+ia+s, wc conclude that 
9^11, 9^S, w and v satisfy the assumptions of the lemma. Since X — S* is not 
incremented, the induction hypothesis applies and proves that there is a constant 
cj with the required polynomial dependence such that 

< C-j |||n - '^tiq,K+lft,P+S + "^7 \\w - v\\n-p,k+i-p,s ■ 
To estimate the second term on the right of (A.7), we fix 

r ^ max{||w||^^^_^i s, \\v\\^j^^^s} , 

set Vj = Br^{0), and recall from above that 11,1] G Cn,k{{Vj},U,I;{Zj}), cf. 
Lemma A. 8. Then, d^w and d^v map U x I into each of the domains Vs, • • . , Vk 
of n, S. Applying the induction hypothesis to 11, S and dfj_w, df^v, we obtain that 
there exists a constant cg such that 

||n o df,W -So a^w||Ar_p_5_i,K-P-S 

< C8 |||n - s|||„_i^^ p^g + C8 ii^^w - a^«IU-i-pK-ps 

< Cg |||n - S|||^_;^^^ Q p^g, + Cg ||w - l'|lAr_p/f+i_ps • 

The first term on the right hand side is estimated as before, yielding a bound of 
the form (A. 9). 

We have thus found the required upper bounds for all terms on the right of 
(A. 6), thereby completing the inductive step also when K is incremented. D 

In the proof of Lemma A. 6, we used part (a) and proved statement (b) of the 
following lemma which we state for later reference. A proof of part (a) can be found 
in [12, Lemma A.7]. 
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Lemma A. 8. Let H, S, w, and v be as in Lemma A. 6; let r > Q, Vj — Br^ (0) for 
j ^S,...,K,0<P < niin(7V - 1, K), and 

Il{w;u, ii) ~ (dwTl o w){u, fi) li and 'S{'w;u, ji) = (dw'^ ° v)(u, (1)111 (A. 8) 
Then: 

(a) II e CAr_i,x.o,s({Vj}, W,I; {Zj}) with |||II|||Ar-i,K,o.s bounded by a polyno- 
mial in \\w\\n-i.k.s o-iT-d r \'^\n.k.o,S; o.'^^d the same holds true for S. 

(b) There is some polynomial c > in |||n|||Ar+i^if_o,s, miN+i.Kfl^s, \\w\\n,k,s, 
\\'>^\\n,k,S7 and r such that 

ll|n - ^\lN-i,K,a,p+s - ^ ll|n - ^\In,k,o,p+s + c\\w - v\\j^_^_pj^_pg . (A. 9) 

Now we are ready to prove the result on the stabihty of fixed points of contraction 
mappings on scales of Banach spaces. 

Theorem A. 9 (Stability of contraction mappings). For N,K ^ Nq with N > K, 
let Z = Zq D 2i D ■ . ■ D Zk be a scale of Banach spaces, each continuously 
embedded in its predecessor, let Wj C Zj be a nested sequence of closures of open 
sets, let X be a Banach space, and let U C X and I C R &e open. Assume 

(i) n,E G Cn+i^k{{W,},U,I;{Wj}); 

(ii) w t-^ n(w; u, 11) and v -^ 5](w; u, fi) are contractions on Wj with contraction 
constant c'- < 1 uniformly for all u £ U , fi €z T, and j — 0, . . . , K . 
Then the following is true. 

(a) The fixed point equation Il(w; u, 11) = w has a unique solution 

weCN+i^K{U,I;{W,}) 

and \\w\\n+i,k is bounded by a function which is a polynomial with non- 
negative coefficients in |||n|||jv+i,x and (1 — c'A^^ . The same holds true for 
the fixed point v — S(w; u, ^). 

(b) Let P < K < N . Then there is some polynomial c with nonnegative co- 
efficients in |||n|||jv+i,x, |||E|||jv+i,x, and (1 — c'j)~^ for j ~ 0,...,K, such 
that 



\\w-v\\j^_pj^_p < c|||n-s|| 



N,Kfi,P ■ 



Part (a) is a version of a contraction mapping theorem on a scale of Banach 
spaces which was already proved in [12, Theorem A. 9]. 

Proof. We prove part (b) only. It follows the same pattern as the proof of part (a) 
with the additional difficulty that we need to carefully keep track of differences in 
the various spaces. As before, we use induction in N and K. For N = K = P, we 
must estimate 

||w-u||2 < \\U{w;u,fi) -n{v;u,fi)\\2 + \\U{v;u,fi) - S(w; u, /^)||^ 

<Co\\w-v\\2 + \\U{v;u,n) -E(u;u,^)||2, 

where Cq is the common contraction parameter with respect to the Zq norm. There- 
fore, 



1 „„ „„ 1 



II"' ^ ^\\ciUxX;Z) ^ T~T l|n ^\\c{WpxUxI:Z) 



lin-sii 



1 _ „/ 11^^ ^\\c(WpxUxi:Z) 1 _ „/ 111^-" ^mp,p,o,p- 

i Cq 1 Cq 
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We first prove that the conclusion also holds when we increment TV, holding K 
fixed. By Lemma A. 3, 

ll^-«llw+i-p.K-p < sup \\{V>uW-V)uv)u\\j^_pj^_p+\\w-v\\^_pj^_p. (A.IO) 
||«||<i 



By induction hypothesis, there is a constant ci which is a polynomial in |||n|||7v+i Ki 



||5^|||Af+i,if J and (1 — c' ) ^ for j = 0, . . . , -ftT such that 



\\w - v\\j^_pj^_p < ci ||n - s|i^^;^^o,p • 

It remains to compute an appropriate bound for the first term on the right of 
(A.IO). 

Note that w{u,^) = D„w(m,^)m, where u = {u,u) & U x B, B = Bi{0), is a 
fixed point of the contraction map 11 given by 

Il{w; (u, u), /i) = dwTl{w{u, ^);u, ^)w + 9„n(w(u, /i); u, /i) u 

= n(u); u, ^) + (9„n(w(u, ii);u,ii)u . (A. 11) 

Similarly, v{u, ii) = D„f(u,/i)u is a fixed point of S. From part (a) we know that 
w,ve Cn+i,k{U,I- {Wj}). Setting 

r = max{|||n|||jv+i,K. 111^1^+1,^} ^J^^^^ Y~d- ^^'"^^^ 

and Vj — Br ' (0) for j = Q, . . . ,K, we find by Lemma A. 6 (a) and Lemma A. 8 (a) 
that n, E G Cm+i,k{{Vj}M X B,I; {Vj}). Hence, IT and S satisfy the assumptions 
of the theorem and, by induction hypothesis, there is some constant C2, depending 
polynomially on (1 — c')^^ for j = 0, . . . , -ftr, |||n|||7v+i,if , and |||E|||7v+i,if such that 

II* - ^\\n-p,k-p < C2 ll|n - S|||^_^o p 

< C2 |||n- E|||jY^gp + C2 \\{duIlow - du'Sov)u\\j^_p^_p 

' (A.13) 

where, in the second inequality, we refer to definition (A. 8) of 11 and S and to 
Remark A.l. By Lemma A. 6 (a) and Lemma A. 8 (a), taking note of Remark A.l, 
the norms 

ll|n|||„+i_^ < \lti\l,^^^K + \\(dunow)u\\^^^j^ 

and |||S|||Ar+i,K are polynomials in |||n|||Ar+2,K, |||I;|||w+2.k, \\w\\n+i,k, \\v\\n+i,k, 
and r. Due to the definition of r in (A. 12) and part (a), these quantities, hence the 
constants in (A.13), have bounds that can be chosen as polynomials in |||n|||7v+2,i<'j 
\imN+2,K, and (1 - c;.)-i for J = 0, . . . , if. 

Applying Lemma A. 8 to the first term on the right hand side of the second line 
of (A.13) and Lemma A. 6 to the second term, both with 5* = 0, we find that there 
is a constant C3 depending polynomially on |||n|||Ar+2,K, |||S|||7v+2,if, II'^^IItv+i k^ and 
II^IItv+i k such that 

II* - ^Wn-P.K-P ^ C3 |||n - mN+j^,K.,0,P + C3 \\w - v\\n-P,K-P 
< C3 |||n - ^\lN+l,Kfi,P + "^4 |||n - ^\In,K,0,P 

< C5 |n — ^\\\N+i,Kfi,p ■ 
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In the second inequality we have used the induction hypothesis so that C4 and C5 
are polynomials in |||n|||Ar+2,_ff , |||S|||Ar_|_2,_fs', and (1 — c'^)^^ for j = 0,...,K. This 
concludes the inductive step in N. 

Second, we prove that the conclusion also holds when we increment K < N , 
holding N fixed. Recall from Lemma A. 4 that 

Ik - v\\n-p,k-p+i < ll«^ - '"\\n-p.k-p+i,i 

we will estimate the three norms on the right hand side separately. For the first 
norm note that a translation of the scale with Zj = Zj^i and the induction hy- 
pothesis show that 

II"' - ""Wn-pm-p+i,! ^ C6 |||n - S|IU,K+1,1,P : (A.15) 

where Cg is a polynomial in |||n|||Ar+i,K+i > |||n|||jv+i.x+i^i, |||E|||Ar_|_i^jf_|_i, and (1 — 
c^)-i for j = 0,...,iC + l. 

For the second term on the right of (A. 14), we apply the induction hypothesis 
on the trivial scale, so that is a constant C7 such that 



|w- wIIat^po < crllln- E|| 



N.P.O,P ■ 



For the third term on the right of (A. 14), we note that w = dfj_w and v — d^v are 
fixed points of the respective contraction maps 11 and S of the form 

n(w; u, /i) = duJ^{w(u, ^);u, ^)w + dfj_Il{w{u, ^); u, /i) 

= n(w; u, ^) + d^Ii{w{u, ^); u, ^) . (A. 16) 

By part (a), v,w & Cn,k+i{1^,I, {-Zj}). Setting 

1 



max{|||n|||^^^^^, |||S|||j^^^_^J -J^^'^kT 



j=0,...,K i - Cj 

and Vj ~ Br ' (0) for 7 = 0, . . . , iiT, we find that, by Lemma A. 6 (a) and Lemma A. 8 
(a), n, S G C7v,k({Vj},Z^,I; {Vj}). Hence, n and S satisfy the assumptions of 
the theorem and, by induction hypothesis, there is some constant cs, depending 
polynomially on (1 — c'-)^^, |||n|||7v./i', and |||S|||7v.k such that 

II* - ""Wn-P-I^K-P < C8 |||n - A\N-l.Kfi,P 

< cs |||n-S|||^_;^^op + C8 \\dfj,now -df,Y.ov\\j^_p_^^_p 

(A. 17) 

where, in the second inequality, we refer to definition (A. 8) of H and S and to 
Remark A.l. By Lemma A. 6 (a) and Lemma A. 8 (a), taking note of Remark A.l, 
the norms 

ll|n||lAr,K < llin||lAr,K + ll^^n o w||^^^ 

and |||S|||Ar,K are polynomials in |||n|||Ar+i,_R-+i, |||S|||w+i,K+i, \\w\\n.k+i, \\v\\n.k+i, 
and r. Due to the definition of r, these quantities, hence the constant cs in (A. 17) 
has a bound that can be chosen as a polynomial in |||n|||Ar+iji-+i, |||S|||Ar+iji--)_i, and 



(l-c')-iforj = 0,...,i^. 
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The first term in the second line of (A. 17) is estimated by Lemma A. 8. We 
obtain 

ll|n - mN-l,K,0,P ^ "^9 ll|n - S|IU,X,0,P + ^^9 11"^ ~ '"\\n-P,K-P 

<cio|||n-E|||jv,K,o,p- (A.18) 

Here cg is a polynomial in |||n|||7v+i,K+i > |||n|||Ar,K, |||i;|||Ar+i^^+i, ||w||Ar^;f+i > 
ll'^^llw-i K ^^'^ W'^^Wn K+ii ^^'^ ^^ have used the induction hypothesis in the last 
inequality, with cio a polynomial in |||n|||Ar+i k+i, |||S|||Ar+i k+i, and (1 — c'^)^^ for 
j = 0,...K. 

For the second term on the right of (A. 17), note that the hypothesis of the 
theorem, with K replaced by K + 1, implies that 

d^n,d^j: eCN,K+i,oA{{Wj},u,^;{Zj}) , 

so that Lemma A. 6 applied with S = 1 yields a constant cn which is a polynomial in 

|||n|||Ar,if+i > \ldf,U\lN^K+i,o,i, Ills III Ar,_R-+i, l|w||w_^+i > ||w||^_^_^i_i, and H^ylU^K+i 
such that 

\\{df,Il)ow- {df,T,)ov\\j^_p_^j^_p 

< cii |||5^n - a^Slll^^^^^ p^^ + cii llw - -flU-p.K+i-p.i 

< cii |||n - s|||^^^^^ „ p + cii ci2 |||n - s|||^^^_^^ -^ p 

< ci3|||n-S|||^^^^iop, (A.19) 

where the second term in the third inequality is due to (A. 15), and ci2 and C13 
depend polynomially on the required quantities. Inserting (A. 18) and (A.19) into 
(A. 17) then concludes the inductive step in K. D 
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